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Abstract 

Statistical mechanics and thermodynamics for ideal fractional exclusion stati- 
stics with mutual statistical interactions is studied systematically. We discuss 
properties of the single-state partition functions and derive the general form of 
the cluster expansion. Assuming a certain scaling of the single-particle partition 
functions, relevant to systems of noninteracting particles with various dispersion 
laws, both in a box and in an external harmonic potential, we derive a unified 
form of the virial expansion. For the case of a symmetric statistics matrix at a 
constant density of states, the thermodynamics is analyzed completely. We solve 
the microscopic problem of multispecies anyons in the lowest Landau level for 
arbitrary values of particle charges and masses (but the same sign of charges). 
Based on this, we derive the equation of state which has the form implied by ex- 
clusion statistics, with the statistics matrix coinciding with the exchange statistics 
matrix of anyons. Relation to one-dimensional integrable models is discussed. 
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1 Introduction 



In recent years it has been appreciated that there exists a nontrivial possibility 
of mutual quantum statistics, or "statistical interaction" between distinguishable 
particles. One example is multispecies exclusion statistics |L], which postulates 
that inserting into the system a particle of species b reduces the number D a 
of single-particle states available for species a by some g a b- The quantities g a b 
constitute the exclusion statistics matrix (later referred to as the statistics matrix) 
and are called mutual statistics parameters for a ^ b. 

Another example is multispecies exchange statistics, or multispecies anyons 0, 
|3|. [I . Their definition is the following: (i) when two particles a are interchanged, 
the two-body wave function acquires a phase factor exp[i7ra a J; (ii) when particle 
a encircles particle b, it acquires a phase factor exp[2i7ro; a fe]. The matrix {a a b} 
will be called exchange statistics matrix. 

In fact, there is a connection between the two models, which was first estab- 
lished for a single species: The thermodynamic quantities for a system of anyons 
(in a box) confined to the lowest Landau level (LLL) of a strong external magnetic 
field H are the same as those of exclusion statistics particles, all having the same 
energy (energy of the LLL) upon identification 5 = a |. This correspondence 
was extended to the multispecies model in the particular case when the electric 
charges of all the species are equal; one then has g a b = a a b @]. 

The relation between the two models was also discussed in the context of an 
algebraic approach to statistics (aimed at searching for an algebra of observables 
for identical particles). The algebraic definition of fractional statistics in one 
dimension |7| applies to anyons in the LLL, where the dynamics is effectively 
one-dimensional |§. On the other hand, this definition suggests a realization of 
ID fractional statistics by systems with inverse square interactions J7| (see also 
@), which has been justified by interpreting the thermodynamics of the Calogero 
and Sutherland models |TD|, [□]] as that of systems of noninteracting particles 
12]. The statistical distribution for ID fractional statistics found in this way ||12| 



turns out to be the same as that derived from Haldane's multiplicity formula for 



exclusion statistics 13 



In addition to the above, (single-species) integrable models of the Calogero- 



Sutherland type were directly related to exclusion statistics [|T4], [15], jib] . It is 
therefore natural to expect that there should exist multispecies generalizations 
of those models that would correspond to multispecies exclusion statistics. 

It is generally believed that excitations in the fractional quantum Hall effect 
(FQHE) can be described either as anyons (by a Berry phase argument |T7|]) or as 
exclusion statistics particles (by a state-counting argument [[I], |18|). In particular, 
when excitations with charges of different magnitudes are present — in multilayer 



systems || and in the hierarchical structure [T!| - the multispecies model is 



applicable. Another problem of interest is that of FQHE quasielectrons and 
quasiholes, again viewed as two species. By counting the states, their off-diagonal 
exclusion statistics parameters were recently shown [EH] to be antisymmetric; on 
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the other hand, the exchange statistics matrix is always symmetric by definition. 
How the two descriptions match in this case, is not clear by now. 

Therefore, developing statistical mechanics for mutual statistics is of interest 
both theoretically and phenomenologically. Thus far, some progress has been 
achieved on the issue |13|, [H], |J, generalizing the results of the better-understood 
single-species case, where the thermodynamics was studied for a system in a box 
|| [23| as well as in a harmonic potential ]15|, |24], ^5 . 



This paper is devoted to systematically studying statistical mechanics of par- 
ticles of arbitrary number of species obeying exclusion statistics and calculating 
their thermodynamic quantities. 

In Sec. |]we introduce the single-state grand canonical partition function and 
prove that it factorizes into a product of partition functions corresponding to sep- 
arate species. We infer a general formula for the coefficients of the expansion of 
its associated single-state thermodynamic potential in powers of the single-state 
Gibbs factors. Those coefficients, multiplied by dimensional factors depending 
on the dispersion law of the particles and the dimension of space, give the clus- 
ter coefficients, which are discussed in Sec. |3|. Assuming a certain temperature 
scaling of the single-particle partition functions — which holds for systems of 
noninteracting particles with various dispersion laws, both in a box and in an ex- 
ternal harmonic potential, — we also derive a unified form of the virial expansion, 
calculating explicitly a few first coefficients. 

Section |4] is devoted to a specific case when the statistics matrix is symmetric 
and in addition the density of single-particle states is constant in energy. This 
case can be analyzed completely, yielding a closed form of the equation of state, 
which may be expanded to give both low-temperature and low-density series. In 
Sec. | we discuss various physical examples of exclusion statistics for systems of 
noninteracting particles both in a box and in an external harmonic potential. 

In Sec. |6] we address the problem of multispecies anyons in the LLL at dif- 
ferent values (but the same sign) of charges and different masses (hence different 
cyclotron frequencies) for different species. In Subsec. |6.1| we solve this problem 
exactly, exposing anyons of different species to harmonic potentials with different 
frequencies such that the resulting level spacing is the same for all the species. In 
Subsec. p.2| we derive the equations of state for multispecies anyons in the LLL 
both in a harmonic potential and in a box. It has the form implied by exclusion 
statistics, with the statistics matrix coinciding with the exchange statistics matrix 
of anyons. Relation to exclusion statistics is also demonstrated by a semiclassical 



argument in Subsec. 6.3 



We conclude with remarks on searching for microscopic integrable models for 
multispecies exclusion statistics. 
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2 Partition functions 



The definition of exclusion statistics is the following Jj]]. If G a is the initial single- 
particle Hilbert space dimension for species a and N a the number of particles of 
that species, then (i) the reduced Hilbert space dimension is 



D a = G a - 9ab(N b - S ab ) , 



(2.1) 



and (ii) the many-particle statistical weights (multiplicities) are determined as 

(D a + jV a -l)! 



w = U 



(2-2) 



V NJ(D a -l)\ 

Applying fl2.1p -( pT2"| ) locally in phase space, so that the total multiplicity is 

W = Y[W {i) } (2.3) 

% 

where % labels groups of single-particle states of nearby energy and is 
given by (gj) with the changes N a -> N®, G a -> G® , and D a -> £>W, makes 
the definition of zdea/ fractional exclusion statistics for many species 1T3L M (see 



also a related discussion of systems with internal degrees of freedom [26 



The above equations make the starting point for constructing statistical me- 
chanics. Assuming that does not depend on a,(] one deduces equations for 
the average occupation numbers n$ = Nf£'/G$ in the thermodynamic limit, 



which can be conveniently written as [13 



w 



(0 



where xty = e^ a and 



w 



n 6 (i 



(0 



(i). 



(') 



(2.4) 



(0 



l + ni l) 



v .,. (2-5) 

It is natural to introduce the single-state grand canonical partition function £W 
which determines the occupation numbers as [13 

9 



?7 



(0 



j; 



('). 



dXa 



It was observed that for a symmetric statistics matrix, g a b 
partition function £W factorizes as 



ne 



(0 

a ) 



(2.6) 
#6a, the 

(2.7) 
However, 



3 This holds for a wide class of practically interesting systems (see below, Sec 
the case of being a dependent could be worked out as well; see, e.g., Eqs. ( |6.16[) -( 6.1S ) 
and the discussion thereafter. 
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where the grand partition function £fp for particles of species a is connected with 
the occupation numbers by [see (|2.5| )1 

& = ■ ( 2 - 8 ) 

The relations (|2.5| )- (|2.8| ) were found in Ref. |4j] by explicitly expanding (for 
the case of two species) the functions £W, an d rvf 1 in powers of the Gibbs 
factors {x^}. We observe here that the above relations can be put on a more 
general ground. Namely, in the generic case of a nonsymmetric g a b, from Eq. ( |2. 1| ) 
we find [we drop the superscript (i) from now on] 

P- = -V% , (2.9) 

OXb Xb 

where 

A*= — + -^-. (2.10) 
w a 1-Wb 

Equations ( |2.6| )— ( [27^ ) then imply that D a b relates the partition functions to the 
occupation numbers: 

E 7l « £> o6 = 6- (2-11) 



Substituting (p[) and (|2~T0| ) herein recovers (pD- This shows that 

are consistent with ( |2.4| ), ( |2.5| ) in the generic case (any dispersion law and a 

nonsymmetric matrix g a b). 

Equations (|2.4|), (|2.71), and fl2.8p implicitly determine £ as a function of the 
Gibbs factors {x a }. In particular, expanding 

ln£= £ fk,..k a x^---x k ^ , (2.12) 

where s is the number of species, one can find the coefficients fki...k s one by one - 
We infer a general formula based on explicit calculations for low orders, 

/fci...fc r o...o = (-l) r ™ — 7 Fki.-.kr ) ( 2 -13) 

lln=l K n 

s — r 

where 

Pl.«l>--> \n=l / 

JV_i,q r _l=l 

Here all the numbers k\, . . . , k r are assumed to be different from zero, and the 
PnS and g n 's summed over in ( 2.14Q are constrained as follows: (i) Pi < P2 < 
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. . . < p r -i] (ii) p n 7^ 9nj (iii) h ma Y n °t be p m = q n and q m = p n ; (iv) at least one 
of the g n 's has to be equal to the number that is lacking in the set {p%, . . . ,p r -i}- 
For example: 

•7~fci = 1 ; 

Tk^hs = klg 21 g 31 + hk 2 g2ig32 + hk^gzi + cycl. perm. , (2.15) 

where "cycl. perm." denotes terms obtained by simultaneous cyclic permutations 
of the subscripts of g's and k's. Hence one gets 

f » ° = s§ (i - kJ r) ("«) 

(cf. @) as well as 

/llO.-.O = — (#12 + ^2l) ; 

/210...0 = -(</i2 + 202i)(l - 25-11 - gw)/2 , 

/1110...0 = #215-31 + #21532 + #23#3i + cycl. perm. . (2.17) 

For two species and a symmetric matrix g ab , the coefficients (|2.13|) reduce to those 
obtained in Ref. j4|. 

One can also evaluate the coefficients of the expansions 

ln£ Q = £ f^. ks x^--.x k /; (2.18) 

k\ . . . hg 

the answer is that is obtained from fl2.13|) by replacing jF fcl fcr with kr , 

which in turn is obtained from (|2.14j) by restricting that p n 7^ a for any n (thereby 
fixing all the p n 's). For example, 

•^fcife = ^1521 , Flxki — kzgvz ; 
■Tfateks = k\g2i9zi + hk 2 g2i932 + fokigmgsi , etc. (2.19) 

For a symmetric matrix g ab , we observe an important property 

k a 

fk!...k s = ^-q - k Jk!...k s , (2.20) 

which will be used in Sec. f| when deriving an equation of state. 

3 Cluster and virial expansions 

The formulas of the previous section enable one to compute the equation of state 
provided the single-particle spectrum is known. Let the latter be = +(^\ 
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where the £W's are common for all the species. By summing (|2.12 ) over single- 
particle states m £ = In S = , one then obtains a cluster expansion 

-/Hl= £ b kl ... ks zt---z k / , z a = e^ (3.1) 

with the cluster coefficients 

b kl ... ks = Z[{KP)e-^ a ^f f ki ks , (3.2) 
where K = J2 s n =i k n and 

Zm = T,e- KW ( 3 - 3 ) 



is the (species-independent) "shifted" single-particle partition function, corre- 
sponding to counting the 

temperature (3 as (cf. |25| 



sponding to counting the single-particle energy from the ground state. One has 
z i(/3) = Eie - ^" = e _/3£ « 0) Z[((3). Assume that Z[(j3) scales with the inverse 



Z W) = §T§ • (3-4) 

We will see below that this scaling is relevant in the thermodynamic limit for 
various physical systems. The cluster coefficients (|3.2|) then become 



The expression for particle numbers N a = z a -J^- InS reads 

N a = £ k a b kl ... ks z^---zf . (3.6) 

fei... k s 

Using ( 2.13 ) to determine the cluster coefficients ( |3.5| ), one then deduces from 



J]) and Q3.6|) the "virial expansion" 



Ni 1 • • • N ks 

-P&= 22 Afc i-^[^/( / 5)lfe 1 +-+fe s -i ' (3.7) 



where the dimensionless "virial coefficients" A kl ___ ks can all be evaluated exactly^, 
even if there is apparently no general expression for them. Up to the third order, 

4 Expressions for virial coefficients in terms of cluster coefficients for two species may be 
found in Ref. M. We give here the simplest formula of that kind for three species: Am = 

2&0116101 + 26oii6no + 26ioi6no - 26m, where b klk2ks = bfcife 2 fc 3 /&ioo & Oio 6 ooi ' as wel1 as tlie 
associated formula relating the cluster coefficient to many-particle partition functions: 6m = 
2-Zooi-Zoio-Zioo — ^100^011 — ^010^101 — ZqoiZuq + 2m. 
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they are: 



Mo. 


.0 = 


i; 




^20. 


.0 = 


-2- 2 -\l-2 9ll ) , 




Alio.. 


.0 = 


2- 1 -\g 12 + g 21 ); 




^30. 


.0 = 


( 4 -l-5 _ 2 . g-2-«) _ 


(4- 5 -3-%i(l- 


^210.. 


.0 = 


-4^(^12 + ^21) (2 


-45-11 -512-521) 






+ (<7i 2 + 2g 21 ) 


(1 - 25n - 512) , 



A1110...0 = -2 • 3 (521531 + 521532 + 523531) 

+ 4~ 1/2 ~ 5 (52l531 + 521532 + 523031 + 01201a) 

+ cycl. perm. . (3.8) 

Note that the mixed virial coefficient of order two vanishes for an antisymmetric 
statistics matrix. 

For a symmetric statistics matrix, the last two expressions become 
A210...0 = -(4- 5 - 3^)512(1 - 2511 - 512) , 

A110...0 = 2(4~ 5 - 3~ 5 ) (512513 + 512523 + 513523) • (3.9) 

If in addition 5 = 0, then the virial coefficients of total order three do not depend 
on the statistics parameters. This case will be discussed in detail in the next 
section. 



4 Thermodynamics for a symmetric statistics 
matrix at a constant density of states 

The case of g a b = gba and 5 = 0, in fact, can be completely analyzed. It is also of 
particular interest since it means a 1/(3 scaling for Z[(/3), i.e., a constant density 
of states in energy. The particle numbers ( |3.6|) read 

N a = Z[(P) £ ka f kl ... ks % ■ ■ ■ **■ . (4.1) 

kl...k s K1 H h Ks 

The statistics matrix being symmetric, using Eq. (|2.20|) yields an important re- 
lation 

N a = Z[{(3) In eaU^ for all 6 , (4-2) 

which enables one to use all the formulas from the previous section provided the 
arguments of the functions £ a and w a , the Gibbs factors {xb}, are replaced by 
the fugacities 

To deduce the equation of state, consider the derivatives 

_djm = ^dz L N 1 

dN a ^dN a z b - 1 • } 
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From ( f4.2| ) we find 

dz b = (M~ l ) ba 
dN a Z[((3) 

dz b 



(4.4) 

where the matrix M ab = ^-ln^ can be expressed, with the use of (|2.8|) and 



9[), as 



and consequently 



M ab = —(D^ab , (4.5) 

z b 



(M~ l ) ba = fD ba . (4.6) 



Inserting this into ( f4.4|) and then the result into ( |4.3| ), we obtain 

gig N a i v iV fe 

<9iV a Z((/5) e^/ z i(« - 1 V^Z;(/5) ' 1 } 

where the symmetry of the statistics matrix was used. Integrating ( [4.7D yields 
finally the equation of state 

- « = jjzbj {?r + ■ <4 - 8) 

Due to the scaling Z[(f3) oc (3 , only zero-temperature terms in the last ex- 
pression depend on the statistics parameters. Equation (f4.8|) can be expanded to 
give both low-temperature and low-density (high-temperature) series. In the for- 
mer case, restricting to the perturbative (behaving as powers of the temperature) 
terms, one obtains 

- Q = E = ^-sZ T 2 + -L J2 9abN a N b , Z = (3Z[ {(3). (4.9) 

b ZZ ab 

This in particular implies that the specific heat is C = ^ti 2 sZqT to all orders in 
T, which is the same as that for a mixture of Fermi gases (one should impose the 
condition g > to do the low-temperature expansion). 

In the opposite limit, for small N a /Z[(/3), one can expand ( |4.8[ ) to obtain 

-^^N a + -U aab -- SM )— + ^—^—j , (4.10) 

where Bk are the Bernoulli numbers (£>2 = h, £>4 = etc., B21+1 = 0). Thus, 
just as Eq. ( |3.8j ) shows, it is only the second order virial coefficients that depend 
on the statistics matrix. Eqs. (|4.8|) - (|4.1CI|) are generalizations of the results of 
Refs. []15, ^4], [23|, [25| to the multispecies case. 
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5 Examples 



In this section we consider particular dispersion laws of particles, mostly dis- 
cussing examples which are related to integrable models. For a single species, 
it has been shown that integrable models with inverse square interactions, the 
Calogero [10| and Sutherland [jllj models, are equivalent thermodynamically to 



systems of noninteracting particles obeying exclusion statistics. This naturally 
suggests that there should exist generalizations of those models corresponding to 
multispecies exclusion statistics. It is in this sense that we refer below to the 
"multispecies Calogero and Sutherland models" (which have yet to be discov- 
ered). The equations of state found below for noninteracting particles obeying 
exclusion statistics apply at the same time to the conjectured multispecies inte- 
grable models, thus providing guidelines for their search. 

5.1 Particles in a box 

Consider a gas of free particles in a Z)-dimensional box of volume V. In the 
quasicontinuous spectrum approximation, when the momentum summation can 
be replaced by integration, for particles with the dispersion law e a (p) = + Ap a 
(a 7^ 0) the shifted single-particle partition function reads [^3], 



1{P1 (2^) D T(l + D/2)(A(3) D /° ' 1 J 

satisfying ( |3.4| ) with 5 = D/a — 1. One has Q = —PV, where P is the pressure, 
and the expansion (|3.7| ) becomes the usual virial expansion for a system in a box: 



PP= E ^ I ...k s p k 1 1 ---P k s s , Pa = N a /V, (5.2) 

fel ■ ■ ■ &s 

where the (dimensional) virial coefficients are 

\ fcxH hk s — 1 

Z[{ 1) 

with A^...^ defined in Sec. |]. 



afci...fc s — A kl ...k a { ryr 7 n\ ) ■ (5-3) 



5.1.1 "Multispecies Sutherland model" (D = 1, a = 2, 6 = -1/2) 

For a single species of particles with a quadratic dispersion law e a (p) = p 2 /2m 
in one dimension, it is the Sutherland model that realizes exclusion statistics. 
For a conjectured multispecies version of that model, the dimensional factor 
in ( |5.3|) relating the dimensional virial coefficients to the dimensionless ones is 



V/Z[(/3) = At, where At = J2izf3/m is the thermal wavelength. The single- 



species expressions agree with those obtained in Refs. [27, ^3 . 
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5.1.2 Chiral fields on a circle (D = 1, a = 1, 8 = 0) 



The model of chiral fields on a circle (a generalization of the model |29| to 
the multispecies case) is a field theoretical realization of exclusion statistics, with 
a symmetric statistics matrix, for particles with linear dispersion (above a gap) 
propagating in only one direction: £ a (p) = + vp, p > 0. Due to the latter 
condition, Z[({3) = V/2tivP, half of the value given by ( |5.1| ). Since 5 = 0, the 
equation of state is given by (|4.8|) ; the low-temperature expansion (|4.9| ) agrees 
with that found in Ref. |28] . 



5.1.3 Nonrelativistic particles in 2 dimensions (D = 2, a = 2, 5 = 0) 

Particles with a quadratic dispersion law in two dimensions and a symmetric 
statistics matrix (a conjectured model of Ref. f3C§, see Sec. |7|) also fall into the 



category considered in Sec. [| so that the equation of state is again ( |4.8| ), with 
V/Z' x {(3) = A|. 

5.2 Particles in a harmonic potential 

We now consider systems of noninteracting nonrelativistic particles in a harmonic 
well of strength w in D dimensions. The single-particle partition function is 

zm = (i - e-^r • (5 ' 4) 

Expanding this to the leading order in f3u results in 

Z'M ^ ^ , (5.5) 



thus respecting the scaling ( |3.4| ) with 5 = D — 1 (with D = 1 as the condition 
for the density of states to be constant). This leads to the equation of state (|3.7| ) 
[and to (|4.8| ) for D = 1 and a symmetric statistics matrix]. 

If one considers the correction terms, of order 0(1) in (|5.5| ), they might lead 
to corrections to the virial coefficients of very high order, N and above. However, 
if the virial expansion converges, these corrections are negligible inside the radius 
of convergence. 

Note that for particles in a harmonic potential, there is no well defined volume 
occupied by the gas; therefore, the equation of state is understood as an equation 
relating the thermodynamic potential Q to the particle numbers, temperature, 
and the harmonic frequency us. Note also that the single-particle dispersion law 
is linear (a = 1) here, because the single-particle energy depends linearly on the 
quantum numbers. 
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5.2.1 "Multispecies Calogero model" (D = 1, a = 1, S = 0) 

For a single species with a linear dispersion law, the Calogero model realizes ex- 
clusion statistics. The "multispecies Calogero model" , which would be equivalent 
to a system of noninteracting particles with a symmetric statistics matrix in a 
one-dimensional harmonic potential |[U| , would again fall into the class considered 
in Sec. HI and be governed by the equation of state (|4.8|), with Z[((3) = 1/ f3ui. 



5.2.2 Harmonic potential in 2 dimensions (D = 2, a = 1, 8 = 1) 

For the same model as in |5.1.3| but involving a harmonic well, i.e., a linear disper- 
sion law, whence Z[(0) = l/(/?u;) 2 , the equation of state is the generic Eq. (|3.7| ) 
with the "virial coefficients" ( |3.8| ). 

Another important example, to be considered in detail in the next section, is 
multispecies anyons in the LLL, both in a harmonic potential and in a box. 



6 Multispecies anyons in the lowest Landau level 

The problem of anyons in the LLL, being of relevance, in particular, to the FQHE, 
was originally solved for a single species in Ref. pi. The multispecies version of 
the problem was addressed in Ref. [|J , with the restriction that all the species have 
the same electric charge and mass, and the equation of state in a box was derived. 
Here we revisit this problem, allowing anyons of different species to have different 
magnitudes (but the same sign) of charges and different masses. We derive the 
equations of state in a box and in a harmonic potential, showing how they fit 
into our general framework, and provide a simple semiclassical picture. 



6.1 Spectrum in a harmonic potential 

The single-particle Hamiltonian in a magnetic field and a harmonic potential is 

H=-—dd-u c (zd-zd) + -^zz (6.1) 
m 2 

(z = x + iy, d = d/dz, 5 = d/dz), and its spectrum is 



E 



U + ^) (wt - uj c ) + (n + ^ (wt + oo c ) , (6.2) 



in 



where u c = —eB/2m (eB < 0), u t = ^uj^ + uj 2 , with £,n = 0,1,2,.... The 
quantum numbers i and n are, respectively, the angular momentum and the 
number of the Landau level. The LLL restriction, n = 0, leads to 

Eg = uj t + , (6.3) 

with vj = u t — uj c . 
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In the multispecies problem, the charges e a and the masses m a being different, 
so are the cyclotron frequencies uo ca . Let the harmonic frequencies u a be different 
as well. The many-body Hamiltonian becomes 



H 



N 



E 



a a (ft ^ a ^ta - 

(Jaj&aj ^ca\^aj^aj z ajOaj) ~r ~ ^aj^aj 
Tfl a Z 



(6.4) 



(j numbers particles of a given species, d a j = d/dz a j). Making an ansatz 



— ^/vexp 



x - m a u ta _ 



(6.5) 



yields the Hamiltonian acting on \& 



N- 



Hn — E dajdaj + (^to — UJ ca)z a jd a j + (k>ta + ^ca) z aj9 a j + ^to 



(6.6) 



We now choose the a; a 's such that = cj t a — u ca be the same for all the species. 
For all the o^'s thus defined to be positive (and to tend to zero, rendering the 
system free, when w — > 0), all the c<J ca 's have to be positive, consequently the 
signs of all the charges have to be the same, such that e a B < 0. Then the 
LLL eigenfunctions of the Hamiltonian ( |6.6j ) satisfying the anyonic interchange 
conditions [HI have the form 



N 



n 



^nn 



(6.7) 



(aj)<(bk) 



sym 



where (aj) < (bk) means: a < b or (a = b and j < k), so that each pair is counted 
only once; and symmetrization is performed over the coordinates of particles of 
the same species only. The energy of the state ( |6.7|) is 



E^ + E 



N a (N a - 1) 



-OLn 



+ -J2N a N b a 



ab 



ab 



w . (6.8) 



The LLL spectrum is obtained by letting the £ a /s run from to 00 with the 
restriction £ a j < i a ,j+i- The partition function is 



Zn^.Ns = exp 



' N a (N a - 1) \ 

2^ 7y a aa + - 2^ N aN b a ab (3w 

■.a 1 1 ab ) . 



n 



a L Utl (1 " C 



(6.9) 
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6.2 Equation of state 

Let us now derive the equations of state for multispecies anyons in the LLL both 
in a (species-dependent, as explained above) harmonic potential and in a box 
(labeling the associated cluster and virial coefficients with superscripts uj and V, 
respectively). 

The cluster expansion for the system in a harmonic potential reads 

ln~= £ &£,.. fcs ^---^, (6.10) 

where the cluster coefficients b% kg can be expressed in terms of the partition 
functions in the standard way, by matching ( |6.10[) to H = Z)fc x ...fc s %kx...k s z \ ' ' ' z s s : 



and the partition functions are given by 

We are interested in going to the thermodynamic limit. In accordance with 
the aforesaid, the harmonic frequencies are u>l = [w + oj ca ) 2 — u^ a , where oj ca = 
—e a B/2m a . The thermodynamic limit is understood as w —>■ 0, so that uj 2 ~ 
2zuu ca — > 0. To the leading order in f3w and u a /u) ca , the cluster coefficients are 

(k± + • ■ ■ + k s )(3w' 



C- — " 

^ki...k s = TL T T u \o^fk±...ks\g ab ->a ah ■ (6-11) 



Appearance of the coefficients fk x ...k s (|2.13|) in this expression indicates that 



exclusion statistics is present in the system. Noting that in the case at hand 
Z[(/3) = l//3tu, we see that Z[(/3) and b% k satisfy (|3.4|) and (|3.5p , respectively, 
with = uj ca and 5 = (constant density of states). We can thus conclude 
that the equation of state for LLL anyons in the harmonic potential considered 
coincides with (|4.8|) , generalizing the result of Ref. f25| for a single species. 



Now we turn to the equation of state for the same system in a box. There 
is a general procedure of deriving the equation of state in a box starting from 
the cluster coefficients in a harmonic potential (originally proposed for anyons 



and then put on general grounds [[3^] ) . Here we follow the lines of Ref. [[32 
Assuming that the external potential varies slowly in space, the relation In 2 = 
(3PV valid in a box is replaced by 

lnS = /3 J d D rP(r) , (6.12) 

where P(r) is the local pressure. The latter is evaluated within a small volume 
where the external potential $ ( r ) acting on the particles of the a-th species is 
approximately constant, using the equation of state for particles in a box with 
the chemical potential \i a replaced by \i a — $ a (r): 

pP(r) = I £ {zie-^) kl ■ ■ ■ (z s e-^) ks . (6.13) 

fei... k s 



Substituting this into ( |6.12j ) and comparing the result to Q6.10 ), one gets 



b l k = y fci ^ s hS>( , . (6.i4) 



14 



For D = 2 and $ a ( r ) = ^a^ar 2 /2, this turns into b kl ks = f3V{kimiUj\ + 
• • • + k s m s ujg)b kl fc 3 /27r. In going to the thermodynamic limit, one has m a uj^ ~ 
2m a wuj ca = 2nzupL a , where p^ a = m a uj ca /T[ is the Landau level degeneracy per 
unit area. Hence, we find the relation between the cluster coefficients of LLL 
anyons in a box and in a harmonic potential (the "thermodynamic limit pre- 
scription"): 

bl... ks = V(k lPL1 + ■■■ + k sPLs )(3wbl... ka • (6.15) 
Substituting ( |6.11| ), one finds 

h v v k ipLi + ••• + k s p Ls k 

b k!...k s - V \- k Jki...k s \g ab ^a ab , (6.16J 

hence, the virial coefficients as in (15.2!) are 



„ _ (k 1 + .- + k a -i)\ v [(^) fc --i]n fe ^ 

which implies the following equation of state for LLL anyons in a box: 

PP = E ln f 1 + ] r /PLa / ) • ( 6 - 18 ) 

V 1 ~ iZb^abpbl Pha) 

This equation also shows that the statistical mechanics of the above system 
of anyons is governed by exclusion statistics. Indeed, starting directly from Hal- 
dane's definition ( |2.2| ) and assuming that all particles of the same species have 
the same energy e a (u ca in the case at hand), one easily derives (see also Ref. ||) 
the equation of state which coincides with ( |S.18[ ) upon changes G a /V — > p^ a and 
Qab ®ab- [For a species-independent p^ a = p^, Eq. ( |6.16|) is of course again a 
special case of (pH)|), with = cu ca , Z[{0) = PlV, 5 = —1.] 

It is possible to give a mean-field interpretation [[H], |[ of Eq. (|6.18 ) in the 



spirit of the Chern-Simons model. In that model, the statistics parameters are 
given by 

a a b = , (6.19) 

where <frb is the Chern-Simons flux, whose proportionality to the charge ensures 
that a a b is symmetric^. Smearing out the fluxes and adding them to the external 
magnetic field changes the Landau density p^a into p£a = PLa — Y^b a abPb- Now, 
the quantity under the logarithm in (|6.18|) is l+p a / pf^ al which is what one gets by 
setting all a a bS to zero and replacing p^ a with p^ a . Equation ( |6.18|) is different 



from the one conjectured previously in Ref. ||21|| , although they coincide for a 
species-independent pl- 

5 Note that if the flux comes from a real solenoid, the corresponding formula is a a i> = 
(e a (f>b + eb4>a)/^, differing essentially by a factor of two from the above |33|| . We thank Alfred 
Goldhaber for elucidating the point. 
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6.3 Semiclassical quantization 



Further insight into the nature of exclusion statistics in the system at hand can 
be gained by considering it semiclassically. The classical equation of motion of a 
particle in a magnetic field and a harmonic potential is, in our notation, 

r = -uj 2 y - 2u c t x e z . (6.20) 

Normal modes, obviously, correspond to circles, with r x e z = —fir (eB < 0, so 
we have chosen positive fi to correspond to clockwise rotation) and r = — fi 2 r. 
Substituting into the above brings up two solutions, 

fi = ±u t - w c , (6.21) 

with u t defined as above [cf. (|6.2|) 1. The upper sign corresponds to the lifted 
degeneracy of Landau levels, the lower one to Landau excitations; opposite signs 
of fi mean that the directions of rotation corresponding to the two modes are op- 
posite. The LLL restriction means that only the former, fi = vj mode is excited, 
so that single-particle orbits are circles, the angular momentum is quantized in 
integer units, and the energy is 

E t = tw . (6.22) 

For a many-body system, the orbits are concentric circles and can be ordered by 
increasing value of angular momentum (or radius). Putting a flux <ft inside the 
orbit of a charge e increases the allowed values of the kinetic angular momentum 
by a = e<f)/2ir. Consequently, the quantization condition now is that the angular 
momentum £ a j of a particle is an integer plus additions from all particles whose 
orbits are inside (cf. p4| ): 

iaj = Lj + E a aALj ~ hk) (6.23) 
bk 



[Eq. ( |6.19| ) has been taken into account], where £ a j are integers and 6 is the step 



function defined here as 8(x) = 1 for x > 0, 6(x) = for x < 0, and #(0) = 1/2. 
The total energy is E = Y, a j Ee aj , and substituting ( |6.22j ) and ( |6.23p herein yields 
the correct result 06.81), save for zero oscillations. 



Eqs. ( |6.22j )- QB.23| ) have recently appeared in Ref. [T2^] as defining exclusion 



statistics (in a somewhat different geometry); basically, they are the Bethe ansatz 
equations. 



7 Concluding remarks 

We have worked out in detail the statistical mechanics and thermodynamics for 
multispecies exclusion statistics starting from Haldane's combinatorial formula. 
As of now, two main groups of models are known to realize exclusion statistics: 
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the two-dimensional LLL anyon model and one-dimensional integrable models. 
As for the first one, we have demonstrated that exchange statistics is equivalent to 
exclusion statistics if the charges of the particles are all of the same sign. The case 
when charges of both signs are present, which would be relevant, in particular, 
to quasielectrons and quasiholes in the FQHE, has yet to be investigated. In 
fact, the spectrum (|6.8|) does not lead to the thermodynamics considered for a 
nonsymmetric matrix g a i, and it is unknown whether the LLL spectrum of anyons 
with different signs of charges does so. 

Concerning integrable models, the most interesting question is as follows: 
What would be the generalization of the Calogero and Sutherland models that 
would realize exclusion statistics in the multispecies case? The virtual coincidence 



of the energy levels of single-species LLL anyons and of the Calogero model pT 



suggests that the energy levels (|6.8|) for the multispecies case might at the same 
time be the energy levels of some integrable model. A similar argument applies 
to the Sutherland model, for which the generalized momenta i would still satisfy 
an equation like ( |6.23| ), by the Bethe ansatz argument, but now with a quadratic 
dependence of the energy of free particles on those momenta. 

One should note that different microscopic spectra may result in the same 
thermodynamics. In Ref. p5| , a quantum mechanical spectrum for a single- 



species two-dimensional system was proposed which leads to an equation of state 
identical to that of two-dimensional exclusion statistics particles with quadratic 
dispersion (see Subsec. |5.1.3|) . The principal difference of our approach from 
that of Ref. is that the spectrum ( |6.8| ) is derived from the rule of the kind 
( |6.23|) for the generalized ( "renormalized" |36|) momenta (a generalization of the 
rule ( |fc>.23 ) can also be given for a two-dimensional system |3D||). In terms of 
these momenta the energy has a form corresponding to a free system [cf. ( |6.22| )1, 
and the effect of statistics thus amounts to changing the bare momenta to their 
"renormalized" values determined by equations of the type ( |6.23| ), in accordance 
with the definition of multispecies exclusion statistics by the Bethe ansatz-like 
equations ||28|| . Along the above lines, the spectrum ( |6.8| ) and its two-dimensional 
generalization may be interpreted as the spectra of the systems discussed in 
Subsecs. |5.2.1| and |5.1.3| , respectively J5(J . 

As was recently discussed, the generalization of the Calogero-Sutherland mo- 
del which keeps the interaction two-body and only allows different couplings 
between particles of different species implies constraints on possible values of 
the coupling parameters already at the level of the ground state of the system 
[ |37jj . One might expect that the multispecies generalizations of the Calogero and 
Sutherland models which would reveal generic mutual exclusion statistics should 
involve three-body interactions, in the spirit of the approach of Ref. which 
exploits analogies with anyons. 
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